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mm- Wave  AlGaN/GaN  HFET's 


1  Introduction 

AlGaN/GaN  HFET’s  demonstrate  considerable  promise  for  advanced  RF  power 
sources  for  communications  and  radar  systems.  Experimental  HFET  devices  have 
produced  RF  output  power  on  the  order  of  10-12  W/mm  of  gate  periphery  at  frequencies 
up  to  X-band.  However,  as  frequency  is  increased  to  the  mm-wave  region  RF 
performance  significantly  degrades.  There  are  a  variety  of  physical  effects  that  currently 
limit  the  high  frequency  performance  of  these  devices.  In  particular,  it  is  proposed  that 
charge  non-confinement  in  the  conducting  channel  of  these  devices  can  produce  space 
charge  and  related  effects  that  limit  the  performance  of  the  HFET’s.  This  project 
investigates  charge  non-confinement  phenomena  and  develops  an  analytic  channel  charge 
model  that  can  be  used  in  a  comprehensive  HFET  model  for  use  in  large-signal  RF  circuit 
simulators. 

2  Drift-Diffusion  Model 

The  properties  of  the  AlGaN/GaN  device  have  been  derived  first  from  a  one¬ 
dimensional  model. 

2.1  Model 

Assuming  the  AlGaN  layer  has  a  fixed  charge  layer,  a  fixed  surface  charge  number 
density  nss0  (induced  polarization  charge)  was  assigned  to  the  interface  between  the 
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AlGaN  and  the  GaN  and  has  been  experimentally  determined:  nss0  =4x10  (cm  ). 

We  start  with  the  differential  equation  for  the  charge  balance  at  equilibrium  in  one 
dimension  (1-D).  Here  the  diffusion  and  the  drift  currents  are  set  equal  in  order  to 
maintain  equilibrium. 

dn(x) 

j total  =  0  =  4  An  n(x)E(x )  +  qDn  — — 

dx 

(1) 

where  we  have  assumed  a  fixed  charge  mobility,  /un  and  a  fixed  charge,  q .  The  electric 
field  E(x)  and  the  charge  density  n(x)  are  both  functions  of  the  distance  into  the 
medium.  In  addition,  we  use  Gauss’  Law  for  electrostatics, 

=  -«„w 
dx  s 

(2) 

and  assuming  the  charge  density  is  differentiable  and  the  electric  field  is  twice 
differentiable,  Gauss’  Law  can  be  written  as 

d  E(x)  _  q  dn(x ) 
dx2  s  dx 


(3) 

Substituting  equation  (3)  into  equation  (1)  results  in  a  Riccati  equation  and  the  boundary 
conditions 


dE(x ) 
dx 


+  - 


2  kT 


E2(x)  =  0 


(4) 


The  constants  are  found  through  a  substitution  of  the  boundary  conditions.  At  the 
interface,  x=0,  the 


£(0)  =  £max  and  n(0)  =  nss 


where  the  maximum  charge  on  the  GaN  side  of  the  interface  between  the  AlGaN/GaN  is 
given  by  nss.  The  maximum  electric  field,  Emax  at  the  interface  can  be  found  from 
Gauss’  law  under  the  assumption  that  the  polarization  charge  at  the  interface  is  the  sheet 
charge  density  introduced  earlier,  «  vv0  =4x10  (cm  ) .  Hence,  from  the  solution  for 
the  electric  field  due  to  a  static  sheet  charge  density,  we  find  the  maximum  electric  field 
at  the  interface,  Emax  =-qnSSQ  /  s2  where  e2  is  the  permittivity  of  the  GaN  material. 

Using  the  appropriate  substitutions,  this  equation  may  be  solved  exactly. 
Consequently,  the  electric  field  distribution,  E(x),  inside  the  GaN  material  is  then  found 
to  be 

^■4i?4e44R_ 

(5) 

The  constants  Ai  and  A2  are  determined  from  the  solution  of  the  Riccatti  equation 


E(x)  = 


\2kTC 


q 


A,  =  1  + 


q 


2  kTC 


A0  =  1- 


q 


2  kTC 


Apply  to  solution  to  get  two  equations  in  two  unknowns: 


4C 


A  A 

(a + a  y 


(6) 


Once  we  get  the  electric  field  as  a  function  of  distance,  E(x),  one  can  find  the  charge 
distribution,  n(x),  from  Gauss’  law 

.  .  s  dE(x) 

n(x)  = - 

q  dx 

where  the  peak  charge  distribution  at  the  GaN  side  of  the  interface  is  given  by  n(0)  =  nss . 
Consequently,  a  simple  derivative  of  the  field  yields  the  charge  density  as  a  function  of 


depth  into  the  material.  Using  these  boundary  conditions,  the  constant,  C  is  given  as 
function  of  the  maximum  electric  field  and  the  charge  at  the  interface, 


C  =  E 


2 

max 


q 

2kT 


(7) 


nss  is  the  injection  level  and  as  injection  increases,  C  gets  smaller.  The  constants  Ai  and 
A2  are  also  a  function  of  nss,  so  the  high  injection  effect  should  be  interesting. 


2.2  Sample  Results 

The  parameters  for  the  simulation  of  the  one  dimensional  model  were  drawn  from 
published  literature.  In  the  following  example,  the  electric  field,  the  charge  number 
density  and  the  total  charge  enclosed  are  given  for  a  variety  of  enclosed  total  charge 
values.  The  first  and  largest  total  charge  enclosed  is  given  by  the  precise  distribution 
which  balances  the  polarization-induced  sheet  charge  density  (red/solid  line).  The  second 
trace  is  roughly  0.5*maximum  charge  density  at  the  interface  for  the  given  (fixed) 
polarization  charge  (blue/small  dash).  The  third  trace  is  roughly  the  0.25*charge  density 
of  the  operating  point  at  the  interface  for  the  given  (fixed)  polarization  charge  (black, 
large  dash). 


Interface  Charge  Density:  1.5*10A19  (  C/cmA3) 
"  "  Interface  Charge  Density:  1.1*10A19  (  C/cmA3) 
—  Interface  Charge  Density:  3 .4*  10A1 8  (  C/cmA3) 


Figure  1:  The  Electric  Field  Intensity  for  (1)  constant  polarization  induced  sheet  charge 
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( nss0  =  4x10  (cm  ) )  and  (2)  a  varying  total  charge  in  the  1  -D  cross-section 


Charge  Number  Density 


Cumulative  Charge  Density  in  the  GaN 


Interface  Charge  Density:  1.5*10A19(  C/cmA3) 
“  “  Interface  Charge  Density:  1.1*  10A1 9  (  C/cmA3) 
—  Interface  Charge  Density:  3 .4*  10A1 8  (  C/cmA3 ) 


Interface  Charge  Density:  1.5*10A19  (  C/cmA3) 

'  "  Interface  Charge  Density:  1.1*10A19  (  C/cmA3) 

—  Interface  Charge  Density:  3.4*10A18  (  C/cmA3) 

(a)  Charge  Number  Density  (b)  Total  Charge 

Figure  2:  The  Charge  Characteristics  for  a  constant  polarization  induced  charge,  but 
varying  maximum  charge  density  at  the  interface 


3  Numerical  Solution  of  the  Two-Dimensional  Problem 

The  solution  of  the  one-dimensional  problem  was  adequate  as  a  first  cut  for  the 
problem  solution  and  we  found  that  its  results  matched  well  with  the  measured  results  in 
the  literature.  However,  the  device  represents  a  three-dimensional  problem;  however,  for 
speed  and  convenience,  a  two-dimensional  solution  was  sought  at  this  point.  This  model 
should  lend  greater  insight  into  the  operation  of  the  device. 

3.1  Model 

The  numerical  solution  was  created  first  using  a  finite  element  model  for  Poisson’s 
equation.  Although  this  will  be  the  preferred  route,  the  finite  element  method  required 
absorbing  boundary  conditions  at  the  free  edges  of  the  device.  At  this  point,  the  finite 
element  model  was  not  cooperating  and  was  “shelved”  in  order  to  construct  a  simpler 
model  using  a  finite  difference  approach.  Figure  3  displays  the  geometry  and  the 
electrical  parameters  of  interest  in  this  problem.  Note  there  are  three  layers:  the  AlGaN 
layer,  the  GaN  layer,  and  the  substrate.  Actual  devices  may  be  constructed  with  a  larger 
number  of  similar  layers  and  the  simulation  code  can  accommodate  these  problems  as 
well.  For  transparent  results,  however,  this  simple  model  was  chosen.  The  contacts  for 
the  drain,  the  source  and  the  gate  are  shown  with  heavy  lines  and  an  associated  voltage. 
The  gate  voltage,  Vq,  is  not  assigned  since  it  will  remain  a  variable  in  the  following 
simulations.  In  each  simulation  event,  two  plots  are  show  for  the  two  dimensional  device: 
a  contour  plot  and  a  surface  plot. 
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Figure  3:  The  AlGaN  Device  Geometry  for  Two-Dimensional  Simulations 


In  each  of  the  following  cases,  the  potential  due  to  the  induced,  fixed  surface  charge  at 
the  interface  is  held  to  the  previously  reported  value:  1.6  x  1014  (C/m2),  the  drain  voltage 
was  held  to  10  volts  and  the  source  was  held  to  common.  The  dimensions  of  the  device 
are  given  in  Figure  3.  The  materials  include  the  AlGaN  layer,  the  GaN  layer  and  a  layer 
of  substrate,  assumed  to  be  sapphire  (the  relative  permittivity  of  sapphire  was  set  to  1 1.6). 
At  the  open  regions  of  the  AlGaN,  the  free  space  region  was  terminated  after  two  levels 
of  spatial  discretization  with  a  radiation-like  boundary  condition.  Hence,  the  field  was 
permitted  to  progress  without  reflection.  This  of  course  could  only  be  approximated. 


Height  4nmj  Height  \nm)  _  Height  (rim) 


3.2  Sample  Results 


Figure  4:  Forward  Biased:  Gate  Voltage  held  to  1  Volt 


Figure  5:  Forward  Biased:  Gate  Voltage  held  to  0  Volts 


Figure  6:  Forward  Biased:  Gate  Voltage  held  to  -1.15  Volt 


Electrostatic  Potential 


Figure  7:  Forward  Biased:  Gate  Voltage  held  to  -2.25  Volts 


This  final  case  was  chosen  such  that  the  channel  has  closed.  One  side  of  the  device 
maintains  a  positive  potential  while  the  other  side  is  negative.  Simulations  reaching 
beyond  this  bias  yield  an  even  stronger  division. 


Electrostatic  Potential 


Figure  8:  Forward  Biased:  Gate  Voltage  held  to  -4  Volts 


4  Particle  in  Cell 

The  next  step  in  this  process  is  to  allow  the  free  charges  to  redistribute  in  the  two- 
dimensional  model.  This  process  will  be  simulated  with  the  particle  in  cell  method,  as 
outlined  in  Figure  9. 


Figure  9:  Calculation  Cycle  for  Particle  in  Cell  [Birdsall,  1991] 
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